Rotating Bose-Einstein Condensates with a large number of Vortices. 



V. Subrahmanyam 
Department of Physics, Indian Institute of Technology, Kanpur-208016, India. 

1 June 2002 

We study a harmonically-confined Bose-Einstein condensate under rotation. Vortex lattice configu- 
rations are investigated through a variational approach. Vortices with more than a unit of angular 
momentum are not stable. We explicitly show that the critical rotational frequency is quite small, 
as was reported in the MIT experiment. The width of the condensate in the xy plane and along z, 
and the vortex density are determined as functions of the rotational frequency. The aspect ratio be- 
comes small for large rotational frequencies, making the condensate an effectively two-dimensional 
system. 
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In recent experiments a large number of vor- 

tices (about 150) have been generated in a Bose-Einstein 
condensate of Na atoms. The vortices form a hexago- 
nal lattice, in a quasi equilibrium, showing a variety of 
properties which differ from a equilibrium system. In 
these experiments, a large angular momentum is im- 
parted to the condensate by rotating it close to the 
quadrupolar frequency. In addition, the condensate is 
subjected to an asymmetric harmonic confining poten- 
tial: LOzI'l'K ~ 20Hz, uJxy/^TT ~ 80Hz, the rotation axis 
being along z in the MIT experiment The high an- 
gular momentum states of the Bose condensate are sim- 
ilar to the quantum Hall systems and type-II supercon- 
ductors. If the rotational frequency Q is exactly same 
as in-plane harmonic frequency ujxy, these states can be 
mapped to the lowest Landau level of quantum Hall sys- 
tems 1). 

The angular momentum that is transfered to the con- 
densate is carried by vortices, each vortex carrying a 
quantum of angular momentum, similar to the vortices 
in a superfluid. The vortices have an effective repul- 
sion, which we shall see later from a variational calcu- 
lation. This implies that a vortex with two quanta of 
angular momentum is unstable against separating into 
two unit vortices. When a large angular momentum is 
deposited on the condensate, the vortices form a regu- 
lar lattice in the xy plane. The unit cell area v of the 
vortex lattice will depend on the rotational frequency, 
the harmonic frequencies, and the two-particle scattering 
length of Na atoms, through the parameters character- 
ing the condensate wavefunction. In the previous work 
of Ho 1^ , the energy of the vortex lattice has been calcu- 
lated through a variational approach within the Gross- 
Pitaevskii paradigm, using a Thomas-Fermi approxima- 
tion. The calculated critical rotational frequency does 
not agree with the MIT experiment [Q . We will carry out 
a variational calculation, with a variational condensate 
wavefunction that is similar to the one used in |^ , that al- 
lows us to estimate the critical rotational frequency which 
is in excellent agreement with the experiment, and the as- 
pect ratio of the condensate (the ration of the rms widths 
in the xy plane and along the z-direction) . 

The condensate wavefunction is determined from min- 
imizing the Gross-Pitaevskii functional [Q 
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where the single-particle Hamiltonians along z and in the 
xy plane are H, = + ^uj^z'^, H^y = -^"^ly + 

2j4^xy{^'^ + 2/^); is the angular momentum along 
z direction along which the condensate is rotated with 
a frequency fl; the two-particle interaction strength is 
g = Anh^asc/M in terms of the s-wave scattering length; 
N is the number of atoms in the condensate, and M is 
the mass of the Na atom. The wavefunction ip is nor- 
malized to unity. The Hamiltonian in the xy plane is 
similar to that of a charged particle in a magnetic field if 
eB/Mc — 2uj,j;y, at uj^y — In this case the energy lev- 
els are Landau levels, with the lowest level degenerate for 
angular momentum quantum number m = 0,1,2.. The 
degeneracy gets lifted for Uxy 7^ ^■ 

In the absence of interactions, for g — 0, the low- 
est energy state is a direct product of harmonic oscilla- 
tor wavefunctions along z direction and in the xy plane, 
with uniform densities (modulo the gaussian variation). 
However, for 17 7^ 0, states with nonuniform densities, 
more specifically with density vanishing at points sig- 
nifying vortices, are favorable. These states will be- 
come stable for fl > flc, where flc is the critical ro- 
tational frequency. We will follow in spirit to es- 
timate the critical rotational frequency. Let us write 
ip(r,z) = <I>m(r)x(^);tiy separating the z direction and 
xy plane, where xi^) = exp(— z^7^/2a^), <i>m(r) — 
Cmr"^ exp{im(j))exp{—r'^ (3'^ /2aly); = ^/KJMuj^ and 

Uxy = ^/fij Mujxy are the oscillator lengths, 7, (3 are vari- 
ational parameters, and Cm,Cz are normalization con- 
stants. The wavefunction carries an angular momentum 



and since l^l 



it also signifies a vortex at the 



origin of strength m. For 17=0, the uniform density state 
(modulo a gaussian variation), viz. m — Q has the lowest 
energy for 7 = l,/3 = 1 in the absence of the interac- 
tions, corresponding to the groundstate of the harmonic 
oscillator. For g 7^ 0, the parameters /3,7 will change 
substantially. The energy of the configuration is (in 
units of the oscillator energy fuvxy) 
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where the dimensionless parameter G = Nasd a-z^f^ 
(using experimental values = 5 x 10^, asc = SOBohr 
radii, — 0.9 x lO^Bohr radii, the interaction strength 
is G ~ 1.1 X 10'*). The coefficient of the first term in the 
above equation, Fi =< > is just the rms width 

of the wavefunction in units of a'^y/ P'^;Ti{m) = to + 1. 
The coefficient of the last term, Ta = 27Taly/ 13'^ J \<^\^, 
is the interaction contribution; r2(TO) — 2to!/to!to!2'". 
Minimizing the energy w.r.t. parameters 7 and (3 for 
the cases of m = (no vortex configuration) and m = 
1 (one vortex of unit strength), we obtain 70 = 71 = 
(a;.M,)2/5(2G)-i/5,/3g = 72(c..,/c.,), /S? = 2(3l For 
the experimental values of ojz/ojxy — 1/4, G — 10^ we 
get 7o ~ 0.08, /3o = 0.16 (in contrast, 7 = = 1, when 
g — ior the no- vortex state). And the energy difference 
between the two configurations 
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will be negative for 51 > f2c where the critical rotational 
frequency is ^d^xy = "^iQ^^xy/^z is about 0.077 for the 
experimental values of G = 10'^,bJz/i^xy — 1/4. This 
compares very well with the MIT experiments Q , where 
vortices have been observed for frequencies as low as 7 
Hz, with the harmonic oscillator frequency of 80 Hz. 

The above discussion establishes that for 17 > J7c vor- 
tex configurations are stable against a uniform density 
state. We examined configurations with larger strength 
of vortex, viz. for to=2,3.. A larger m will lower the in- 
teraction term, which is clear from the coefficient, but 
will increase the rms width of the wavefunction. For 
TO = 2, the energy is minimized for 72 = 70(8/9)^/''^, and 
131 = 2/3261/5^ implying K2 - = l.2b-il{{%/Sfl'- - 
\)ujz/oJxy > 0. It can be seen that for to > 2, the energy 
is not favorable, implying a vortex with angular momen- 
tum larger than one unit is not stable. The to = 2 state 
is unstable against a state with two m — 1 vortices sepa- 
rated by some distance. The energy is lower for a larger 
separation, implying an effective repulsion between the 
vortices. We can investigate states with a large angu- 
lar momentum by looking at configurations with a large 
number of vortices with unit strength. This is tanta- 
mount to separating the vortices in state $„. Instead 
of a m-vortex sitting at the origin, we will have m one- 
vortices spread all over. The effective repulsion between 
the vortices gives a finite spacing. A regular lattice struc- 
ture is easily probed, either triangular or rectangular (the 
only Bravais lattices in two dimensions). 

Let us consider a trial wavefunction [p| , with q vortices 
each carrying a unit of angular momentum 
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where we introduce a vortex function f{z) = Yii^ ~ 
zi),zi = xi + iyi is the complex coordinate of the I'th 
vortex (carrying a unit of angular momentum) . We have 
also written the wavefunction as a superposition of var- 
ious configurations $,„ (carrying one vortex of strength 
to). Noting that J ^mHxy^m' = 0, for m 7^ to', the 
energy of this state can be written as before (Eq.2) with 



and 
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The evaluation of the above coefficients Fi, F2, is quite in- 
volved, as the computation of the integrals is not straight- 
forward, because of the presence of the vortex function 
f{z). We shall try to extract the leading order behaviour 
of these coefficients as functions of the ratio TTa^y/ j3'^v, 
where v is the unit cell area of the vortex lattice. 

Firstly, we note that the vortex function can be written 
as f{z) = z{z^ - 6f)(z6 - b^){z^ - 6|)(zi2 - 6f ).. , for a 
triangular lattice (actually a hexagonal lattice, including 
the z direction) where bn is the nth neighbor distance in 
units of ttxy/P- Each bn occurs with a power which is 
equal to the coordination number of the n'th neighbor 
vortex. For triangular lattice the coordination is either 
6 or 12. In the case of a square lattice the coordination 
number is either 4 or 8. Here we have fixed the origin on 
one of the vortices. In the product {z — zi)..{z — zq) of 
the first neighbor vortices, only two terms are nonzero, 
viz. z^ , ziZ2..ze = —bf. Similarly for the other neighbour 
vortices. We now resort to a simple approximation to 
estimate f{z). Let us write 



log/(z) = ^log(z-z,) 



R <.27r 

r'dr' / d9log{z~r'e''^), (6) 
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where we replaced the sum by an integral, and 1^ is the 
unit cell volume of the vortex lattice. In the above ap- 
proximation, namely replacing the density of vortices by 
a constant, we have ignored the finite-size corrections 
that arise, as the lattice size is not too big (about 150 
vortices or less). This approximation is similar to the 
average- vortex approximation of ||^ . The details of finite- 
size effects on the vortex function will be given else where. 
Here, we will keep only the smooth part of the density of 
vortices. In the integral the upper limit, viz. the size of 
the vortex lattice, is taken to be R, with ttR^ /v = q, the 
total number of vortices. The integral is easily done, by 
doing an integration by parts in r' variable, keeping the 
log function as the second function, we get (using \z\ = r) 



logf^^qn + \og{^r^^ + ^. 
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Now, = {qv/irY exp {■kt'^ /v), and in turn the con- 

densate density is given by 
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In the limit of ^ iOxy, 0^ tends to one, and 7 will also 
approach unity, as can be seen from the above equation. 
However, even for ^/ojxy = 0.9, the right hand side of the 
equation is quite large (because G ~ lO'^, ujxy/uJz — 4), 
implying 7 is quite small. Here, we can approximate 
4j- — 7^ w (this is tantamount to dropping the kinetic 
energy term in Hz (Eq.l), as usually done in the Thomas- 
Fermi approximation [^), and thus obtain (valid for fl > 



1 - a = 



nr'^ /v 



/32(2G)i ^ to. 



Figure 1: The log of the vortex function is plotted against 
r^, both from exact numerical computation and from 
Eq.7. The size of the lattice is about R=6.3, contain- 
ing about ttR^/u ~ 145 vortex lattice points. 
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In the above a is the ratio of the area of the effective Lar- 
mor circle and the unit cell area of the vortex lattice. In 
Fig.l we plot the numerically computed function log 
for a a lattice with 145 sites, along with the approximate 
value that we calculated by doing an integral. As can 
be seen from the figure, the above computation (Eq.7) is 
in excellent agreement with the numerical computation, 
implying the finite-size effects are quite small here. 

The total condensate density jV'P has now three vari- 
ational parameters, namely, a the number of vortices in 
an effective Larmor circle, (3 the inverse rms width in the 
xy plane, 7 the inverse rms width along z direction. The 
coefficients ri,F2 are calculated as 



1 



As n 



^xy, 



^~ (2G)t(l-/34)t 

—^ 1 , a — > 1 , which would imply 7^1, 
which implies the rms width along z direction will tend 
to a^. In the neighborhood of < tOxy, 7 has the form 
shown in the above equation. Similarly, as — > 0, = 
/3g = 0.16,7 = 70 = 0.08, a = 0. The rms widths of the 
condensate in the xy plane, and along z direction are 



Rxy = V< r2 



> = 



*'xy 



< > = 



V27' 



In the quantum-Hall limit (17 = LOxy), the unit cell area of 
the vortex lattice becomes ly = Tra^y, implying as many 
vortices as the number of particles in the condensate. 
The width of the condensate becomes very large in the 
xy plane, and the width along z direction becomes very 
small (and saturates at a^/v^), as can be seen from the 
above relations. 
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In the limit of a — > 0,ri = r2 = 1, we have the no 
vortex state. Now, our task is to minimize the Gross- 
Pitaevskii energy (Eq.2) over variation of the parameters 
a,/9, 7. Setting the derivatives of i^(a,/3, 7) w.r.t. the 
three variational parameters to zero, we obtain 
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(^-7') = G7/3^(l-a) 
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Table 1: For various values of the rotational frequency, 
the number of vortices within a size of the rms width of 
the condensate density, the rms width of the condensate 
in the xy plane, and the rms width along z direction. 
We have used the experimental values of G = 10*, and 
i^z/i^xy = 1/4. Note that within a size of l.bRxy (the 
condensate density falls to 10 per cent) the number of 
vortices is 2.25 times that shown in the second column; 
for example the number of vortices in the condensate 
is effectively about 25 for the second case shown. As 
n increases the rms width along z direction decreases, 
which will saturate at az/\/2 in the limit Q ujxy 



The first two equations imply = fl/ujxy, and the last 
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The aspect ratio Rz/Rxy becomes very small in this 
Hmit, and the system is essentiaUy two- dimensional, in 
analogy with the quantum-Hall system. The excitations 
of the system are in xy plane. For f2 = ujxy we should go 
back to the original problem of interacting Bose particles, 
rather than working within a Gross-Pitaevskii paradigm, 
as small perturbations hitherto neglected will become im- 
portant. 

The optimal density of vortices can be written as 
1 a/32 2Am 



which reduces to the familiar l/i/ ~ 2Mil/h in the limit 
of SI ^ LUxy The optimal number of vortices within a 
size of the rms width of the condensate is 

TT < > a 
V 1 — a 

The optimal Gross-Pitaevskii energy is given by 

5 2s 7 , UJ, ^1 , a^r ^ 2 , ^ 2 

K = -_iV= ^ 5 ^ 5 1 _ —y.hoJxy. 

4 7r5 uJxy ujiy 

In Table 1 we have listed the optimal values of the varia- 
tional parameters, the rms widths, the number of vortices 
within the rms width, for various values of the rotational 
frequencies. These results are different from the work 
of |3|. For Q^jiOxy ~ 0.9, there are a large number of 
vortices, and 7 is still small, the density profile along z 
direction is similar to the inverted parabola of , except 
near the edge where it has a gaussian tail. But, as f2 be- 
comes closer to LOxy^ as 7 — > 1, the density profile along 
z direction is that of the harmonic oscillator groundstate 
(the rms width is equal to az). A quantitative measure 
of how far the system is from quantum-Hall limit can be 
taken as 1 — 7. 

In the foregoing analysis, the leading order behaviour 
is valid for ujxyi where a large number of vortices 
are generated. The lattice structure of the vortices en- 
ters only through the density. The triangular lattice is 
favored over the square lattice, as for the same density 
of vortices {i' = \^c^^/2 = c^^, where ctr,Csq are the 
lattice spacing for the triangular and square lattices re- 
spectively), the spacing is larger in the case of triangular 
lattice. Since the vortices have an effective repulsion, a 
larger spacing is favored for energy minimization. But, 



in the quantum-Hall regime {fl = LOxy), the energies of 
the two lattice structures become degenerate (square lat- 
tice energy approaches that of the triangular lattice from 
above). Then it becomes important to consider all per- 
turbations, though small and hence hitherto neglected, 
to lift the degeneracy to find the true groundstate con- 
figuration. This implies one should go beyond the Gross- 
Pitaevskii for the condensate wavefunction. There is also 
a possibility of formation of higher angular momentum 
vortices (for > ujxy), and maybe an instability of frag- 
menting of the condensate. This is because, once the 
unit cell area becomes equal to the area of Larmor circle, 
the angular momentum that can be carried by m = 1 
vortices gets saturated. If a larger angular momentum 
is transfered to the system, m > 1 get generated, or the 
system breaks into smaller pieces. It will be interesting to 
investigate other variational states in this regime, as the 
physics here is similar to that of a Quantum Hall sys- 
tem. We are investigating these issues currently, along 
with the question of multi-component Bose condensates 
where the internal degree of freedom introduces many 
features motivating many novel variational states. 

In conclusion, we have calculated the critical rotational 
frequency ilc/i^xy 0.07, above which vortices become 
stable, which compares very well with the MIT exper- 
iments. We have shown that variational states, with a 
gaussian density profile along z direction, and density 
profile with vortices of unit strength , the locations form- 
ing a triangular lattice, in the xy plane get stabilized for 
larger SI. We have computed the width of the wavefunc- 
tions along z direction, and in the xy plane, the aspect 
ration becoming very small as f2 ^ uJxy The unit cell 
area of the vortex lattice becomes equal to the Larmor 
circle area in this limit, implying as many vortices as the 
number of condensate atoms. 
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